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Abstract 

We show that, when boosted to the infinite momentum frame, the quark and gluon orbital angular mo¬ 
mentum operators defined in the nucleon spin sum rule of X. S. Chen et al. are the same as those derived 
from generalized transverse momentum distributions. This completes the connection between the infinite 
momentum limit of each term in that sum rule and experimentally measurable observables. We also show 
that these orbital angular momentum operators can be defined locally, and discuss the strategies of calculat¬ 
ing them in lattice QCD. 
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Apportioning the spin of the nucleon among its constituents of quarks and gluons is one of 
the most challenging issues in QCD. The quark spin measured in deep-inelastic scattering (DIS) 
experiments contributes about 25% to the proton spin [[ll], and a lot of experimental and the¬ 
oretical effort has been devoted to determining the rest pieces in the past 25 years. A recent 
global analysis [0] that includes the high-statistics 2009 STAR JH] and PHENIX [0] data shows 
evidence of non-zero gluon helicity in the proton. At = 10 GeV^, the polarized gluon dis¬ 
tribution Ag{x, Q"^) is found to be positive and away from zero in the momentum fraction range 
0.05 < X < 0.2. However, the result presented in [Q] still has large uncertainty in the small-x 
region. Given that the integral value of Ag{x, Q'^) from x = 0.05 to 0.2 is about 40% of the proton 
spin [ 0 ], there is still room for substantial contribution from the quark and gluon orbital angular 
momenta (0AM). To fully understand the proton spin structure, one needs to find the spin and 
0AM operators whose matrix elements give the partonic contributions that are both measurable 
in experiments and calculable in theory, especially, in lattice QCD. In this work, we address the 
question of orbital contributions, and show that the 0AM operators in the gauge-invariant sum 
rule for nucleon spin of Chen et al. [Bi , when boosted to the infinite momentum frame (IMF), 
are equal to those derived from generalized transverse momentum distributions (GTMD) that can 
be measured in hard exclusive processes. The operators in B 0] can be defined locally on the 
lattice, and therefore our work paves the way for a well-defined non-perturbative calculation of 
the measurable quark and gluon 0AM. 

The quark spin and gluon helicity measured in DIS experiments can be incorporated into the 
naive proton spin sum rule by Jaffe and Manohar B]: 

- = 2 ^^ + + Cg , (1) 


where each individual term is defined to be the proton matrix element of canonical spin and 0AM 
operators in the IMF (or on the light-cone) with the light-cone gauge. As is well known, AS and 
AG are the first moments of polarized quark and gluon distributions, which are given by light- 
cone correlation functions and consistent with operator-product expansion (OPE) [SI] . Meanwhile, 
there are attempts to define canonical 0AM distribution functions and relate them to experiments. 
One recent proposal based on quark models suggests that Cg is related to a transverse momen¬ 
tum distribution (TMD) function hj-j, iffiil] which is measurable in semi-inclusive DIS experi¬ 
ments [Il2l.ll3n. More recently, in a model-independent proposal, Cg is directly given by a Wigner 
distribution or GTMD Ff 4 ifldl - fl^ : 


Cq{x) 


j X ki_ywlc{x,h^,ki_) 

- j d^k^-^Fl^x, 0, kl, 0, 0) , 


( 2 ) 


where M is the nucleon mass, b± and k± are the relative average transverse position and mo¬ 
mentum of the quarks, and “EC” stands for a specific choice of the gauge link that corresponds 
to semi-inclusive DIS or Drell-Yan processes. There is also similar relationship between Cg and 
the gluon Wigner distribution or GTMD. The definition of Cg in Eq. dH) has clear part onic inter¬ 
pretation, and in principle it should be measurable in hard exclusive processes [I19I. I20I1, although 
currently no proper process is identified and the existence of Fi 4 has been doubted [ 12 ill . 

Since the OAM’s from GTMD are given by non-local operators on the light-cone, they are 
not directly accessible in lattice QCD. Nevertheless, there has been recent progress on the lattice 
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calculation of TMD functions ||22| - |24|1 , where a spatial gauge link is used in the calculation and the 
result is evolved to the IMF to obtain the light-eone correlator. A similar formalism is proposed 
in the framework of large-momentum effeetive theory (LaMET) |2^27], where GTMD—along 
with other parton physies—can also be ealeulated in lattice QCD with this approaeh [l25n . This 
can be one direction of calculating Cq and Cg, and one has to faee the ehallenge of renormalizing 
non-local operators with Wilson lines on a finite lattice, whieh was attempted in |2^ 2^, but is 
not yet eonelusive. 

A different perspeetive is provided in as the gluon spin operator in [Bi is shown to be 
equal to that defined from the first moment of ^g{x, Q'^) in the IMF limit, thus it ean be used for 
the lattiee ealeulation of AG. This idea is further developed in [|29ll . and an attempt to ealeulate 
AG in lattiee QCD has been earried out very reeently j^ lsill. 

Inspired by this diseovery, we explore the eonneetion between GTMD and the 0AM operators 
in the gauge-invariant deeomposition of nueleon spin by Chen et al. 


J = J X —i/j + j d^x'ip^x X iDgavsiJ 

+ / d^x E X Aphys + / d^X E\x X Dpure)^ 


phys 5 


(3) 


where i is the spatial Lorentz index. Here and ] are the 

gauge-eovariant derivatives aeting on the fundamental and adjoint representations, respectively, 
with A^ = -f A[^yjg. To make each term in Eq. ([3]) gauge invariant, one requires that under a 
gauge transformation g{x). 


3lphys(a^) ^ 9{x)A^i^y^{x)g {x) , 


^pure(a^) ^ 9{x)[ A^„,,(a;) --9^)9 


(4) 


In addition, to find a solution, it is suggested la that Apj^^^ satisfies the non-Abelian Coulomb 
eondition, 

= , (5) 

while A^yjg is eonstrained by a null-field-strength eondition, 

+ = 0. (6) 


According to [l28h . in the IMF limit A+^g = A+, so one obtains a first-order linear equation for 
A^yj.g from Eq. dS]), 

a+Ap^’^ - ^7/“'‘=A+’'>Ap^’^g = . (7) 

The solution is given by 




( 8 ) 


where the light-eone eoordinates = {x^ ± x^)l\/2. The kernel /C(^' — ^ ) is related to the 
boundary eondition Q. Here £ is a light-eone gauge link defined in the adjoint representation. 
With A^j^“ , we ean easily obtain 


Ap^^;(r) = A'^’“(r) + T^ / dr/c(r-r)5M+’^(r)^^“(r,r). 


(9) 
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After integration by parts, this solution is identieal to the one found by Hatta [16, 321. 


) = - 7 ^ I dy -X , 


xy- 


'yx 


( 10 ) 


where the light-eone gauge link is defined in the fundamental representation. Note that it 
only requires = A+ or = 0 in the IMF limit to get Eq. (|7]), so Eq. (jS]) actually belongs 
to a universality class of conditions that can be used to fix and [{2^. 

Under an infinite Eorentz boost along the z direction, the equal-time plane is tilted to the light- 
cone, and the bilinear operator ^1 • • • ^ transforms into '4)'^^ • • • The quark 0AM Lq defined in 
Eq. dl]) becomes 


lim L'= / 


P^^OO 


pure / 


( 11 ) 


where the covariant derivative on the light-cone with A[^^“ (,^“) given by Eq. (|8]). 

In the same limit, with z = 1, 2 transforms into F^~^, = F^~ remains the same, 

whereas Apj^^^ becomes A+j^^^ that vanishes according to 112811 . Therefore, the IME limit of the gluon 
0AM Lg in Eq. dU) is 


Now recall the quark GTMD function, 

fix,hr, At) ^ I 

^)\PS) , 


( 12 ) 


(13) 


zy ^ is a straight-line gauge link 


where P = P - A/2, P' = P + A/2, A+ = 0, and _ 

2 ’ 2 

connecting the transverse fields at light-cone infinity. Here ±(X) correspond to the link choices for 
semi-inclusive DIS and Drell-Yan processes. 

InlfTbll. it is proved that 

‘'dx£kTky(x,kT,AT) = iE{p's\my(P^-i'5y)m\ps), (M) 


where with A^^^^g given exactly by the solution in Eq. dS]). Therefore, 


pure 


lim {LA = 


P^^OO 


2P+ (2ir)33(3)(0) 


{ps\ [drd^(Tmi*p’e(Ay-E^,)m\ps) 


= lim 

A ^>0 


d 


/ dx(fkT /c^/(x, fcr, At) 


d'^k±-^Flr^{x, 0, fci, 0, 0) , 


(15) 


with = 1, and = 0. 
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Similarly, we can prove that 


lim {LD 

P^^OO ^ 


2P+ (27r)35(3)(0) 


{PS\ I drd^^TFliOe^^CV^^ur.A-^M)\PS) 


= e*-^ lim 


d 


A-i^o idAL 


dxd^kx kj.g{x, kx. At) , 


( 16 ) 


with the gluon GTMD [[I^,[33ll. 


g{x,kT,AT) = 


2xP+ 


dz d‘^ZT i^p+^-_ 
(27r)3 


ikr-ZT ^p/ 




±=o,in(Y.y)l^’S), (17) 


where the gauge links are defined in the adjoint representation. 

This finishes our proof that the quark and gluon 0AM in Eq. ([3]) in the IMF limit are the same as 
the gauge-invariant definitions from GTMD functions. Furthermore, when fixed to the light-cone 
gauge, they reduce to the canonical 0AM in the Jaffe-Manohar sum rule. 

Since the operators in Eq. ([3]) are independent of real time, they can be calculated on the Eu¬ 
clidean lattice. With the EaMET developed in [|^ . one can relate their lattice matrix elements 
to parton spin and 0AM in the IMF through a perturbative matching condition, which has al¬ 
ready been derived at one-loop order Odil . As far as lattice calculation is concerned, solutions for 
Apjjyj, and M satisfying Eqs. ([4}® can be obtained through a gauge link fixed in the Coulomb 
gauge OSL 13611. One starts from the gauge link U^{x) = exp{—iagoA'^{x)) that connects x to 
X + afi. Under a gauge transformation g{x). 


U^{x) —)■ U'^{x) = g{x)U^{x)g ^{x + afi) . 


(18) 


By finding a gauge transformation gc that makes U^{x) = gc{x)UZ{x)g^ ^{x + a(x) where UZ{x) 
is fixed in the Coulomb gauge, one can define a new gauge link and obtain the solution for 


<hys m. 


Upure = 9 c{x)gP{x + ap), 

^phys = — (U^ix) - (^p^ure(^)) = —9c{x){U^{x) - l)c/"^(x) + 0{a) 

= gc{x)Ac{x)gZ\x) + 0(a) . (19) 

One can check that so defined satisfies the gauge transformation law in Eq. (jH) with being 
unchanged and gc transforming as g'^ = ggc- A short calculation confirms that Eqs. dSHb]) are also 
satisfied up to 0(a) corrections which vanish in the continuum limit, 

^pure(^) = ((^p^ure(a^)(^pure(a^ + + ok)UF^^(x) - l) -f 0(a) 

a go 

= 0(a) , (20) 

^^*^phys(a^) = -^9c(x) {Ul(x) - Ul(x - ai)\gZ^(x) + 0(a) 

Qj \ / 

= 0(a) . (21) 
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In this way, is constructed locally on the lattice. The scheme has been applied in the 
recent attempt to calculate AG It can also be used to calculate the 0AM in Eq. ([3]). 


Note that Lq resembles the mechanical 0AM in the sum rule by Ji Il37h . except that the co¬ 


variant derivative in the latter is replaced by Ji’® based on a gauge-invariant 

and frame-independent decomposition of the nucleon spin, so each individual contribution is the 
same in arbitrary frame including the IMF. This allows for a lattice calculation of the quark 0AM 
and total gluon angular momentum in the rest frame, which has been carried out on a quenched 
lattice OSll . In practice, the explicit space coordinate poses a problem for direct calculation of the 
forward matrix element of 0AM due to periodic boundary conditions ll^ . so in ll^ the quark 
0AM is not directly calculated. Instead, it is obtained by subtracting the quark spin from the 
total quark angular momentum. The latter is calculated from the Belinfante-improved symmetric 
energy-momentum tensor , 


Jl = 2 ^*'' j > 


( 22 ) 

(23) 


The angular momentum Jq is obtained from the forward limit of two form factors of Tq'', 


j, = ^\Tm+uo)] . 


(24) 


One might think of the same approach to obtain by using the symmetric energy-momentum 


tensor 


t;""- = i 


(25) 


However, with this definition, Eq. (1221) cannot be satisfied; that is, the right-hand side of Eq. 
will produce the quark angular momentum plus extra terms. Therefore, the sum rule in Eq. 
cannot be satisfied. If, however, one adopts the asymmetric energy-momentum tensor. 


rrasv = 


pure 


Y'lp , 


(26) 


one will obtain Lq using Eq. (1221) . The parametrization of the off-forward matrix elements of a 
similar asymmetric tensor that uses instead of has been discussed in where there 
are two additional form factors, which is not difficult for lattice calculations. However, the frame- 
dependence of spoils the simple parametrization in lldOll . so one has to introduce a temporal 
vector = (1, 0, 0, 0) to include more Eorentz structures and form factors, which is analogous 
to that formulated on the light-cone 114ill . If all the form factors can be calculated in lattice QCD, 
then the quark 0AM will be obtained through a relevant sum rule. 

Alternatively, one can directly calculate the 0AM from the off-forward matrix element by 
utilizing the relation 




{P'S\ 


d^x^p^xWl^^^^p\PS) = 


Umei,{P'S| j , (27) 


and the challenge here is to have a reliable extrapolation to the —)■ 0 limit. 

The third lattice approach is the calculation of Cq and Cg from GTMD which can be generalized 
from the lattice study of TMD as mentioned in the introduction [|22l-l24ll. One of the challenges 
here is the renormalization of non-local gauge-link operators on the lattice. 
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In conclusion, we have proved that, in the IMF limit, the gauge-invariant quark and gluon 0AM 
defined by Chen et al. JH, are the same as those from GTMD. The former ean be ealeulated 
through loeal operators on the Euelidean lattiee, whieh reduees the difficulty caused by non-loeal 
Wilson lines in GTMD, and their matrix elements are matehed to the physieal quark and gluon 
0AM through LaMET. With the development of GTMD measurement in hard-exelusive proeesses, 
we will eventually be able to eompare the theoretieal and experimental results on the parton 0AM. 


This work is partially supported by U.S. Department of Energy grants DE-SC0013065 and 
DE-FG02-93ER-40762. One of the authors Y. Z. is espeeially thankful for the hospitality of the 
Physies Department of the University of Kentueky during his visit when the work is made possible. 
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